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INTRODUCTION 


A  number  of  computer  programs  for  analyzing  effectiveness 
of  fragmenting  warheads  are  currently  in  use.  Many  of  these 
require  a  description  of  the  fragmentation  pattern  and  dis¬ 
tribution  of  velocities  of  fragments  by  angular  zone,  as  measured 
with  respect  to  the  warhead.  The  ra w  data  upon  which  this 
information  is  to  be  based  can  be  obtained  by  firing  the  warhead 
in  a  test  arena  consisting  of  halves  of  two  (nearly)  cylindrical 
surfaces,  one  made  up  of  (flat)  "distribution  panels"  to  gather 
fragment  density  data  and  one  made  up  of  "velocity  panels"  used 
for  recording  fragment  velocity  data.  Because  all  fragments  are 
assumed  to  have  equal  weights,  the  data  consists  of  observable 
fragment  holes  in  the  (distribution)  panels  as  well  as  films  of 
fragments  penetrating  through  the  (velocity)  panels  of  the  arena. 

The  computer  program  ZONE  has  been  developed  in  order  to  provide 
a  convenient  means  of  converting  this  raw  data  into  the  required 
quantitative  description  of  the  density  and  velocity  distributions 
of  the  fragments  by  angular  zone. 

Use  of  ZONE  and  the  methods  described  in  this  report  enable* 
one  to  conveniently  record  the  results  of  the  test  firing  on 
data  cards  without  the  need  of  recovering  any  fragments,  and  it 
provides  a  systematic  approach  to  obtaining  a  description  of  the 
warhead's  fragmentation  pattern  by  angular  zone.  This  can  be 
easily  done  for  arbitrary  arrangements  of  the  zones  without  the 
need  of  any  fragment  data  modification  whatsoever,  and  with  only 
a  very  slight  change  in  program  input.  In  contrast  to  other 
programs  of  its  type,  ZONE  allows  one  to  collect  the  data  using 
rectangular  panels  but  processes  it  to  obtain  an  accurate  picture 
of  the  fragment  pattern  by  conical  zone  through  its  use  of  steradian 
measure. 

The  first  two  sections  of  this  report  describe  the  test  arena 
(and  actual  test)  and  the  method  of  calculation  used  by  ZONE.  The 
description  of  the  arena  is  a  general  one  and  is  not  intended 
to  be  highly  detailed  because  the  methods  of  setting  up  the 
arena  and  obtaining  the  data  are  basically  standard;  more  compre¬ 
hensive  descriptions  of  these  procedures  can  be  found  elsewhere. 

The  final  two  sections  of  the  report  describe  the  computer 
program  ZONE  (including  an  input  guide  and  output  guide)  and  a 
sample  case. 
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1.0  DESCRIPTION  OF  THE  TEST  ARENA 


A  top  view  of  the  test  arena  is  illustrated  in  figure  1. 
The  dimensions  given  in  the  figures  are  those  used  in  the 
accompanying  sample  case  {see  Sec.  2.0);  neither  the  test  nor 
the  computer  program  require  that  these  some  dimensions  be 
used. 


The  warhead  is  placed  in  the  center  of  the  arena,  which 
is  also  the  center  of  the  two  semi-circular  cylindrical  surfaces, 
one  consisting  of  "distribution  panels"  (also  called  "density 
panels")  and  one  consisting  of  "velocity  panels"  (for  obtaining 
velocity  information) .  These  will  now  be  described  individually. 

1.1  Distribution  Panels 

The  distribution  panels  are  usually  made  of  rectangular 
pieces  of  1/2  inch  plywood  on  which  a  grid  is  drawn  dividing  the 
panels  into  sub-panels  or  targets.  After  the  warhead  has  been 
fired,  the  fragment  holes  in  each  sub-panel  are  counted  and 
recorded.  The  horizontal  indexing  system  for  identifying  the 
sub-panels  is  indicated  in  fig.  1  (also  in  fig.  2)  which  shows 
the  distribution  panels  as  viewed  from  inside  the  arena.  The 
letters  H  (Horizontal)  and  V  (Vertical)  are  used  to  label  the 
sub-panels;  for  example,  the  sub-panel  containing  the  asterisk 
in  fig.  2  would  be  referred  to  as  sub-panel  V-7,  H-4  or  as  the 
( 7 , 4 ) th  sub-panel. 


Once  the  number  of  panels  is  decided,  one  can  choosi  a 
radius  so  that  these  panels  will  very  nearly  fit  along  a  semi¬ 
circle  with  this  radius.  The  "effective  radius"  R  of  the  arena 
is  then  computed  as 

NW  Half  circumference  of  arena  (1) 


where  N  is  the  number  of  distribuiton  panels  horizontally,  and 
W  is  the  width  of  each  panel. 1  That  is,  since  the  panel  width 
is  small  compared  to  the  radius  of  the  arena,  one  can  approximate 
the  circumference  of  the  (half)  arena  by  NW.  The  two  ends  of 
each  panel  should  be  placed  at  this  distance  R  from  the  burst 
point  (measured  on  the  ground) ;  this  also  ensures  proper  angular 
placement  of  the  panels.  Although  the  plywood  panels  are  flat 
rather  than  curved  and  the  ends  of  the  sub-panels  will  not 
generally  lie  exactly  on  the  circumference  of  the  semi-circle 
(unless  the  sub-panels  and  panels  are  identical) ,  they  will 
adequately  approximate  a  semi-cijole  of  radius  R. 


^  The  program  actually  uses  NSWS  instead  of  NW,  where 
Ns  and  Ws  are  the  corresponding  quantities  for  sub¬ 
panels;  however,  it  is  obvious  that  NSWS=NW. 

2 

The  actual  radius  of  a  semi-circle  on  which  the  ends 
of  each  panel  would  lie  exactly  is  given  by  (cont.  «~  ) 
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1.2  Velocity  Panels 

The  velocity  panels  are  made  of  aluminum  (with  black  backing) 
attached  around  the  edges  to  a  plywood  frame.  These  panels  are 
also  considered  to  be  divided  into  a  grid  of  sub-panels  which 
are  numbered  as  in  fig.  1.  Just  before  the  warhead  is  detonated, 
a  number  of  long  duration  flash  bulbs  are  set  off  and  the  cameras 
pick  up  the  light  coming  through  the  holes  made  by  the  fragments 
as  they  penetrate  the  aluminum.  Each  hit  is  read  off  the  film 
and  the  location  and  frame  number  of  the  hit  are  noted.  The 
location  is  again  described  by  the  horizontal  and  vertical  indices 
of  the  sub-panels,  which  are  now  formed  by  putting  a  grid  over 
the  developed  film  rather  than  by  drawing  them  on  the  panels 
as  was  done  for  the  distribution  panels.  This  data  together 
with  the  camera  speeds  and  drag  information  allows  one  to 
compute  the  average  velocity  of  fragments. 

The  radius  of  the  velocity  panel  half  of  the  arena  is 
computed  in  a  fashion  analagous  to  that  of  the  distribution 
panel  half. 

It  is  critically  important  that  the  warhead  lie  parallel  to 
and  above  the  center  line  of  the  arena  and  that  it  have  zero 
degree  azimuth  and  elevation  angles.  This  is  because  the  zones 
are  defined  relative  to  the-  direction  the  shell  is  pointing,  which 
is  assumed  to  be  parallel  to  the  center  line  of  the  arena. 


(cont . ) 


R  = 


W 


2 


sin 


/  O  ' 
\  '-/ 


The  difference  in  results  from  using  'vl)  in  place  of  (2)  is 
generally  negligible  (and  for  reasons  of  practicality  it  is 
usually  sufficient  for  the  arena  to  be  set  up  with  a  radius 
within  two  inches  of  the  value  calculated  using  (1)  ) . 
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2.0  METHOD  OF  CALCULATION 

The  purpose  of  ZONE  is  to  obtain  estimates  of  the  number 
and  velocity  of  fragments  by  angular  zone,  from  the  input  data 
gathered  at  the  test  arena.  The  program  is  divided  into  two 
major  sections,  one  for  density  calculations  and  one  for  velocity 
calculations.  Each  of  these  will  be  considered  separately. 

2 . 1  Density  Calculations 

Aside  from  the  basic  dimensions  of  the  arena,  the  data 
necessary  for  the  computation  of  densities  are  the  number  of 
fragment  hits  counted  for  each  distribution  sub-panel.  These 
densities  to  be  calculated  are  called  "fragment  densities"; 
they  are  measured  in  units  of  number  of  fragments  per  steradian 
and  they  provide  a  description  of  the  distribution  of  the  fragments 
from  the  warhead  throughout  the  arena .  They  should  not  be  confused 
with  the  "density  of  a  fragment"  which  might  be  measured  in 
grains  per  cubic  centimeter  and  is  a  characteristic  of  the 
fragment  rather  than  of  the  fragmenting  warhead. 

The  program,  assuming  that  each  panel  is  divided  into  an 
integral  number  of  sub-panels  and  that  these  lie  on  the  perimeter 
of  the  arena,  computes  the  coordinates  of  the  four  corners  and 
the  center  of  each  sub-panel.  Using  these,  it  then  calculates 
the  steradians  (area  on  the  unit  sphere)  of  each  target  or  sub¬ 
panel  and  which  zone  it  is  in.  The  angles  defining  the  zones 
are  inputs  to  the  program,  and  may  be  changed  from  run  to  run. 

The  program  repeats  the  above  for  each  sub-panel,  keeping 
track  of  the  total  fragments  in  each  zone  and  total  steradians 
of  the  targets  in  the  zones.  After  going  through  all  the  sub¬ 
panel  data  the  program  computes  fragment  density  (fragments/ 
steradian)  for  each  zone,  and  under  the  assumption  that  the 
fragment  distribution  is  uniform  wi thin  any  one  zone,  multiplies 
this  by  the  total  steradians  in  the  zone  to  obtain  the  total 
predicted  number  of  fragments  for  each  zone. 
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A  zone  is  defined  as  the  region  of  space  between  two 
cones  having  the  same  axis  and  vertex,  described  by  an 
upper  and  a  lower  angle.  A  point  is  in  a  particular 
zone  if  the  angle  between  the  vector  from  the  burst  point 
to  that  point  and  the  vector  from  the  burst  point  through 
the  nose  of  the  shell  lies  between  the  upper  and  lo''er 
angles  of  the  zone.  The  program  assigns  all  the  hits 
in  any  one  sub-panel  to  .‘hat  zone  which  contains  the 
center  of  the  sub-panol.  This  method  greatly  simplifies 
the  calculations  without  significant  loss  of  accuracy. 
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2.2  Velocity  Calculations 

The  purpose  of  these  calculations  is  to  convert  the  fragment 
data  from  the  velocity  sub-panels  into  a  breakup  of  average 
velocity  by  zone.  As  mentioned  previously,  the  data  is  obtained 
from  high  speed  films  of  the  actual  test,  from  which  frame  numbers 
corresponding  to  each  hit  are  recorded  and  are  inputted  to  the 
program  along  with  camera  speeds  (in  frames  per  second) .  Since 
there  will  usually  be  four  cameras  filming  the  test,  four  camera 
speeds  must  be  inputted  together  with  the  number  of  horizontal 
sub-panels  each  one  filmed. 

As  was  the  case  for  density  calculations,  the  data  is 
collected  by  sub-panel;  here  however,  for  input  purposes  sub¬ 
panels  with  no  hits  should  be  omitted.  The  program  computes 
average  velocities  for  each  fragment  by  using  the  following 
relations: 

Frame  Number  (3) 

T=  Camera  Speed 

D  (4) 

v-  T 

where  T  =  time  required  for  fragments  to  hit  the  velocity  panels 
D  =  distance  from  the  burst  point  to  the  center  of  the  sub-panel 
V  =  average  velocity  of  a  fragment  in  travelling  distance  D. 

The  initial  velocity  of  a  fragment  at  the  burst  point 
immediately  after  detonation  is  computed  using  the  formula 


The  quantity  of  (in  units  of  ft”-*-)  is  a  drag  factor  which  is 
a  slowly  varying  function  of  velocity  and  a! so  depends  on  the 
fragment  shape  and  weight.  The  formula  for  alpha  is 

ol  ?=  CD  j>  k  m"i//3  (6) 

where  C'c  drag  coefficient  (dimensionless) 

?  =  density  of  air  (lb/ft^) 

k  =  shape  factor  (fi;2-grains  i//3/lb) 
m  =  fragment  weight  (grains) 
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P  (RHO)  is  assumed  to  have  a  value  of  .0765  (density  of 
air  at  15 °C)  and  is  included  in  a  DATA  statement  in  the  program. 
Values  of  m  (GRAINS) ,  Cp  (CD)  and  k  (SHAPE)  are  inputs  to  be 
chosen  by  the  user.  A  set  of  typical  values  is  m  =  30.0  ,  Cp  = 
.55  and  k  =  .462,  which  yields  a  value  of  of  =  .006239.  Values 
for  CD  and  k  can  be  found  in  reference  1.  Capitalized  names  in 
parentheses  above  are  the  FORTRAN  names  of  the  variables  preceding 
them.  In  the  program,  &  is  called  ALPH  and  the  term  in  square 
brackets  in  (5)  is  called  FCTR . 

The  program  stores  the  number  of  fragments,  the  sums  of 
the  average  ana  initial  velocities,  and  the  sum  of  the  squares 
of  the  initial  velocities  for  each  zone,  from  which  it  computes 
an  "average  average"  velocity  (AVERAGE  VBAR) ,  an  average  inital 
velocity,  and  a  standard  deviation  (for  initial  velocities  only) 
for  each  zone.  These  are  calculated  as: 

AVERAGE  VEAR  FOR  Ith  zone  =  1  V1  (7) 

nTTT  Z  VIJ 

AVERAGE  INIT.  VEL.  for  Ith  zone  =  1  V  (8) 

nTTTZ.  V0Tj 

STD  DEV  for  Ith  zone  =  j  1 

|n(i)-i 

where  _ 

V  =  average  velocity  for  the  Jth  fragment  in  the 

IJ  Ith  zone, 

V°  =  initial  velocity  for  the  Jth  fragment  in  the 

IJ  Ith  zone. 

N(I)  =  number  of  fragments  in  the  Ith  zone. 

All  summations  are  on  J,  from  J  =  1  to  J  =  N(I). 
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3.0  COMPUTER  PROGRAM 


ZONE  is  a  computer  program  written  in  FORTRAN  for  use 
on  CDC  6000  series  computers;  it  can  easily  be  modified  for 
use  on  other  computers.  The  next  sections  contain  some  notes 
on  the  program,  descriptions  of  the  input  and  output,  and  a 
source  listing. 

3 . 1  Program  Notes 


1.  Present  limits  in  the  program  are  35  zones  (36  zone 
angles),  16  sub-panels  vertically,  and  90  sub-panels  horizontally. 
These  can  be  increased  by  varying  dimensions  at  the  beginning 

of  the  main  routine  and  possibly  varying  some  output  formats. 

2.  The  inputs  NDENS  and  NVEL  indicate  to  the  program 
which  calculations  are  to  be  performed,  and  the  variable  NLEFT 
will  adjust  the  results  ; ..  correspond  to  the  direction  the  warhead 
is  facing  (fig.  1  corresponds  to  NLEFT  =  1) .  This  does  not 
affect  the  numbering  of  the  sub-panels  ,  but  it  does  reverse 

the  positions  of  the  zones,  since  they  are  measured  relative 
to  the  direction  the  warhead  is  fa  ing. 

3.  The  coordinate  system  used  to  locate  tne  corners  of 
the  sub-panels  is  left-handed  with  origin  at  the  burst  point, 
y-axis  perpendicular  to  the  ground  and  pointing  up,  and  x-axis 
along  the  warhead  axis  pointing  towards  the  far  edge  of  the 
last  panel. 

4.  The  velocity  section  cf  the  program  resembles  the 
density  section  in  that  some  quantities  occur  in  both  sections 
which  may  have  different  values  in  the  two  sections.  In  these 
cases  the  suffix  "2"  is  appended  to  their  FORTRAN  names  in  the 
velocity  section.  For  example,  whereas  R  represents  the  radius 
of  the  density  half  of  the  arena,  R2  represents  the  radius  of 
the  velocity  half. 

5.  If  one  is  mainly  interested  in  zones  totalling  less 
than  180°,  say  from  60°  to  120°,  it  is  advisable  to  include 
two  extra  zones,  from  0°  to  60°,  and  from  120°  to  180°,  to 
account  for  the  remainder  of  the  fragments.  If  this  is  not 
done,  fragments  not  in  the  60°  -  120°  range  will  not  be  printed 
in  the  summary  and  no  count  or  record  of  them  will  be  outputted. 

6.  The  purpose  of  SUBROUTINE  PROJ  is  to  compute  the  area 
of  a  sub-panel  in  steradians  (using  the  coordinates  of  its 
four  corners)  and  to  compute  what  zone  the  sub-panel  is  in 
(using  the  coordinates  of  its  center) . 
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3. 2  Input  Guide 


The  following  is  a  list  of  input  cards  necessary  to  run 
ZONE.  The  input  data  is  divided  into  three  groups  of  cards  I, 

II  and  III,  and  card  sets  in  these  groups  are  A,  B,  C.... 

The  user  has  the  option  of  selecting  density  and/or  velocity 
calculations  as  indicated  to  the  program  by  the  values  of  NDENS 
and  NVEL.  These  values  govern  which  groups  of  cards  are  necessary 
for  the  run.  A  table  showing  which  inputs  are  required  for 
the  different  calculations  follows  the  description  of  the  input 
cards . 

Contents,  Description  Format 

I .  GENERAL  INPUTS 

A.  (TITLE (M),  M=l,  3)  (8A10) 

TITLE  is  any  (80  character  max.)  title  to  be  printed 
at  the  top  of  each  page  of  output. 

B.  NDENS, NVEL, NLEFT, NSET, BRSTH  (4I3,F8.2) 

{0  Do  not  calculate  densities 
1  Calculate  densities 

{0  Do  not  calculate  velocities 
1  Calculate  velocities 

{0  If  warhead  is  facing  last  panels 

1  If  warhead  is  facing  first  panels  (as  in  fig.  1) 

NSET  =  Number  of  zones 
BRSTH  =  Burst  height  (ft) 

C.  (C0NA (K) ,  K=1 ,  NSETl)  (8F10.3) 

These  are  the  cone  angles  (in  degrees)  describing  the  zones. 
The  zones  are  assumed  to  be  contiguous,  so  that  the  number 
of  angles  required  will  be  one  more  than  the  number  of 
zones . 

II.  DENSITY  CALCULATION  INPUTS 

A.  BASE,  WTHSP,  HTSP  (4F10.1) 

BASE  =  distance  above  ground  of  bottom  of  density  panels  (ft) 
WTHSP  =  width  of  density  sub-panels  (ft) 

HTSP  =■  height  of  density  sub-panels  (ft; 
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B.  NV,NH 


(213) 


These  are  the  numbers  of  density  sub-panels,  vertically 
and  horizontally  respectively. 


C.  ( (N(I,J) ,  1=1, NV) ,J=1,NH)  (2613) 


N(I,J)  is  the  number  of  fragments  in  the  (I,J)th  sub-panel, 
read  in  one  (vertical)  column  at  a  time.  Each  new  column 
must  begin  on  a  new  card. 


III.  VELOCITY  CALCULATION  INPUTS 

A.  BASE2,  WTHSP2,HTSP2, GRAINS,  CD,  SHAPE  (6F10.1) 

The  first  three  of  these  are  the  same  as  in  II-A,  except 
they  are  for  velocity  sub-panels. 

GRAINS  =  weight  of  fragment  (grains) 

CD  =  drag  coefficient  (dimensionless) 

SHAPE  =  shape  factor  (ft2  -  grains1''  /lb) 


B.  NV2,  NH2 


(213) 


These  are  the  same  as  in  II-B,  except  they  are  for  velocity 
sub- panels . 

C.  (FSD(M),  KL (M) ) ,  M=l,4)  (4 (F10 .1 , 15) ) 

FSD(M)  =  camera  speed  (frames/sec)  of  camera  M 
KL(M)  =  number  of  sub-panels  (horizontally)  filmed  by 
camera  M,. 


D.  I,  J,  NN,  (KF (M) ,  M=1,NN)  (313,1X1415/(1615)) 

I  =  index  of  sub-panel,  vertical  direction  (V) 

J  =  index  of  sub-panel,  horizontal  direction  (H) 

NN  =  number  of  fragments  in  (I,J/th  sub-panel. 

KF (M)  =  frame  number  of  Mth  fragment  in  (I,J)th  sub-panel. 

Ordering  of  fragments  within  a  sub-panel  is  arbitrary,  as  is 
the  ordering  of  the  input  cards  in  III-D.  A  sub-panel  with 
no  fragment  hits  should  not  be  inputted.  If  there  are  more 
than  14  frame  numbers  to  be  included  for  any  sub-panel,  the 
additional  frame  numbers  must  be  included  in  succeeding  cards 
using  (1615)  format. 
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The  program  may  be  run  for  the  density  half  of  the  arena, 
the  velocity  half,  or  for  both  halves.  The  input  variables 
NDENS  and  NVEL  in  I-A  indicate  to  the  program  which  calculations 
are  desired.  The  table  below  indicates  the  inputs  necessary  in 
the  three  cases: 


NDENS 

7 

0 

1 


NVEL  Necessary  Input  Data 

0  I,  II 

1  I,  III 

1  I,  II,  III 


Type  I  data  is  always  required  and  the  input  cards  are 
ordered  as  they  appear  each  letter  (e.g.,  III-B,  III-C) 
indicating  at  least  one  new  card. 

3.3  Input  Data  Forms 

Convenient  forms  for  reporting  the  data  to  the  program  user 
have  been  prepared  and  are  illustrated  m  the  following  three 
pages.  It  is  a  simple  matter  to  punch  the  input  data  in  the 
correct  format  once  the  data  has  been  entered  on  these  forms 
correctly. 

Form  1  contains  the  title  or  description  of  the  round, 
the  burst  height,  and  the  necessary  test  arena  data  for  the  dis¬ 
tribution  and  velocity  sub-panels.  It  also  contains  the  camera 
speeds  of  the  four  cameras  and  the  corresponding  indices  of  the 
sub-panels  they  film.  A  "REMARKS"  section  has  been  included 
for  any  other  pertinent  information  about  the  test. 

Form  2  is  used  to  record  the  number  of  fragment  hits  for 
each  distribution  sub-panel,  and  as  indicated  above,  H  represents 
the  horizontal  index  of  the  sub-panel  and  V  the  vertical  index. 
The  H  numbering  goes  from  left  to  right  and  the  V  goes  from 
bottom  to  top  .<='•'  that  the  panels  appear  on  this  form  as  they 
are  in  the  arena.  This  data  corresponds  to  data  set  II-C  above. 

Form  3  is  used  to  record  the  film  data  from  the  velocity 
panels.  The  vertical  index  of  the  sub-panel  is  entered  in  the 
first  column  (labelled  V),  the  horizontal  index  is  entered  in 
the  second  column  (labelled  H) ,  and  the  number  of  fragments 
hitting  a  particular  sub-panel  is  entered  in  the  third  column. 

The  frame  numbers  of  each  fragment  arc  recorded  in  the  re¬ 
maining  space.  If  there  is  not  enougn  space  on  one  line  for 
all  the  frame  numbers,  succeeding  lines  can  be  used.  Sub-panels 
with  no  hits  should  not  be  recorded.  The  data  on  Form  3 
corresponds  to  data  set  III-D. 
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FORM  NO.  1 


Title  or  Description  of  Round 
Burst  Height 


Density  Data  Included?  YES  NO  (Circle  One) 

Dimensions  of  Distribution  Sub-Panels:  ,  Width 


Number  of  Distribution  Sub-Panels: 

Distance  Above  Ground  of  Bottom 
of  Panel 

Velocity  Data  Included?  YES 

Dimensions  of  Velocity  Sub-Panels: 

Number  of  ’^lo^ity  Sub-Panels: 


Height _ 

Horizontally 

Vertically 


NO  (Circle  One) 

Width _ 

Height _ 

Horizontally _ 

Vertically _ 


Distance  Above  Ground  of  Bottom 
of  Panel 


Camera  Speed  (Frames/Sec) 

No.  of  Horizontal  Indices  Covered 

Camera  1 

Camera  2 

Camera  3 

Camera  4 

REMARKS : 

11 


FORM  NO.  2  -  DISTRIBUTION  SUB-PANELS 


CO 


I  “A 


FORM  NO.  3  -  VELOCITY  SUB-PANELS 


7^--.  —  '4.--V * ■• ._ '^y.i « ?•«<  ^yjr^V&Tf.r*?^'  fm > l-'^rl-SSsr 1 ll 1 1  lJ. 


3.4  Output  Guide 

The  first  page  of  the  density  output  contains  a  listing  of 
the  dimensions  of  the  distribution  sub-panels,  including  all  the 
information  about  these  found  on  Form  1.  It  also  contains  a 
list  of  the  inputted  zone  angles. 

The  next  page(s)  contains  a  listing  of  the  number  of  fragments 
by  sub-panel,  40  lines  per  page,  and  corresponds  to  the  data 
entered  on  form  2. 

The  final  density  page  contains  a  summary  of  the  results  by 
zone,  including  definitions  of  zones  by  angles,  fragments  counted, 
density  of  fragments  per  steradian,  and  predicted  number  of  frag¬ 
ments  computed  using  the  method  of  Section  2.1.  The  "TOTAL 
NUMBER  FRAGMENTS  PREDICTED"  at  the  bottom  of  this  page  is  simply 
the  sum  of  the  predicted  numbers  for  each  zone. 

The  first  page  of  the  velocity  output  contains  a  listing  of 
the  dimensions  of  the  velocity  sub-panels  (including  all  the  infor¬ 
mation  about  them  found  on  form  1) ,  the  drag  information,  the 
camera  speeds  and  the  inputted  zone  angles. 

The  next  page(s)  contains  a  listing  of  the  number  of  fragments 
by  sub-panel  (40  lines  per  page)  which  corresponds  to  the  data 
entered  in  the  "V",  "H"  and  "NO.  HITS"  columns  on  Form  3. 

The  final  density  output  page  contains  a  summary  of  the  results 
by  zone,  including  definitions  of  zones  by  angles  and  fragments 
counted  per  zone.  The  column  labelled  "AVERAGE  VBAR"  is  an  "aver¬ 
age  average"  velocity  for  each  zone,  since  the  velocities  being 
averaged  are  already  average  velocities  over  the  distance  from 
the  burst  point  to  the  velocity  panels.  The  "AVERAGE  INIT.  VEL. " 
and  the  "STD.  DEV.  INIT.  VEL."  for  each  zone  are  computed  as 
explained  in  Section  2.2. 

If  the  program  reads  any  frame  numbers  equal  to  zero,  it 
will  print  out  the  vertical  and  horizontal  indices  of  the  sub¬ 
panel  and  the  number  of  the  fragment.  For  example,  if  the 
message 


1=3  J=G  M=7 

appeared,  it  would  indicate  that  the  7th  frame  number  for  sub¬ 
panel  V-3,  H-6  was  equal  to  zero.  This  usually  means  that  the 
number  of  fragments  for  this  sub-oanel  punched  on  the  input  card 
(type  III-D)  was  greater  than  it  should  have  been.  This  message 
allows  one  to  locate  certain  incorrectly  punched  input  cards. 
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4.0  SAMPLE  CASE 


The  arena  used  in  the  sample  case  is  shown  in  fig.  1; 
it  is  set  up  so  that  the  panels  are  centered  about  a  height  of 
10  ft.,  corresponding  to  the  burst  height.  The  distribution 
panels  are  16  ft.  high  and  4  ft.  wide,  and  the  bottom  panels 
are  2  ft.  off  the  ground.  There  are  27  of  these  panels, 
each  divided  into  sixteen  2'  x  2'  sub-panels,  so  that  there 
are  54  sub-panels  horizontally  and  8  sub-panels  vertically. 
Using  equation  (1)  the  radius  of  this  half  of  the  arena  is 
R  =  34.3775. 

The  velocity  panels  are  8  ft.  high  and  4  ft.  wide,  and 
the  bottom  panels  are  6  ft.  off  the  ground.  There  are  17 
of  these  panels  each  divided  into  thirty-two  1'  x  1’  sub¬ 
panels,  so  that  there  are  68  sub-panels  horizontally  and  8  sub¬ 
panels  vertically.  Using  equation  (1)  the  radius  of  this  half 
of  the  arena  is  R  =  21.6451. 

For  the  sample  case,  the  values  of  m,  Cq,  and  k  mentioned 
in  2.2  are  used,  resulting  in  the  value  cxf  =  .006239. 

Cameras  1,  3,  and  4  filmed  four  panels  each  and  camera 
2  filmed  five  panels.  This  is  reflected  in  the  entries  on  form 
no.  1  which  lists  the  number  of  sub-panels  horizontally  that 
each  camera  filmed. 

The  sample  case  is  further  described  by  the  data  on 
completed  forms  1,  2  and  3  which  follov’  in  4.1;  these  are 
simply  a  tabulation  of  the  fragment  hits  and  camera  results 
from  an  arena  as  illustrated  in  fig.  1.  A  list  of  the  input 
deck,  card  by  card,  is  given  in  4.2,  and  the  output  from  this 
same  sample  case  is  provided  in  4.3. 
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4 . 1  Completed  Input  Forms  for  Sample  Case 


FORM  NO.  1 


Title  or  Description  of  Round 
Burst  Height 

Density  Data  Included?  i^YESj  NO 

Dimensions  of  Distribution  Sub-Panels: 

Number  of  Distribution  Sub-Panels: 

Distance  Above  Ground  of  Bottom 
of  Panel 

Velocity  Data  Included?  £ye1T)  NO 

Dimensions  of  Velocity  Sub-Panels: 

Number  of  Velocity  Sub-Panels: 

Distance  Above  Ground  of  Bottom 
of  Panel 


CiV-C 


1  o.o ' 

Width 

(Circle 

...  ✓ 

&  , 1 

One) 

Height 

2- O' 

Horizontally 

S' 4 

Vertically 

> 

8 

2. o' 

Width 

(Circle 

hO  ' 

One) 

Height 

1.0 ' 

r 

Horizontally  * 

68 

Vertically 

8 

6.  o' 

REMARKS:  m  -  30  ?ra'*s 

VJorKCftcl  Mt  Ca: 
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SUB-PANELS 


i 

i 

V 

H 

NO. 

HITS 

FRAME  NUMBERS  OF  HITS 
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7 

IS 
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1 

73  ;  . 
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26 
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2-7 

3 

6  3,  J\9S 

FORM  NO.  3  -  VELOCITY  SUB- PANELS 


FORM  NO.  3  -  VELOCITY  SUB-PANELS 


FORM  NO.  3  -  VELOCITY  SUB-PANELS 


4.2  Input  Cards  for  Sample  Case 


The  following  pages  contain  a  listing  of  the  actual  data 
cards  used  for  the  sample  case.  The  numbers  in  the  far  left 
column  are  card  numbers;  they  are  listed  here  for  easier 
reference  and  should  not  be  considered  part  of  the  actual  data 
punched  on  the  cards  which  begins  to  the  right  of  the  line. 
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Output,  from  Sample  Case 
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